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We present a measurement of the branching fraction and time-dependent CP asymmetry of 
B° -> p°K°. The results are obtained from a data sample of 227 x 10 6 T(4S) -> BB decays 
collected with the BABAR detector at the PEP-II asymmetric-energy B Factory at SLAC. From a 
time- dependent maximum likelihood fit yielding 111 ± 19 signal events we find B(B° — > p°K°) = 
(4.9 ±0.8 ±0.9) x 10" 6 , where the first error is statistical and the second systematic. We report the 
measurement of the CP parameters S p o K o = 0.20 ± 0.52 ± 0.24 and C p o K o = 0.64 ± 0.41 ± 0.20. 



PACS numbers: 13.25.Hw, 12.15.Hh, 11.30.Er 

Decays of B° mesons to the p°K° final state are 
expected to be dominated by b — > s penguin ampli- 
tudes. Neglecting Cabibbo-Kobayashi-Maskawa (CKM) 
suppressed amplitudes, the mixing-induced CP violation 
parameter S p o^o should equal sin 2/3, which is well mea- 
sured in B° — > J/ipK° decays 0. Within the Standard 
Model (SM), only small deviations from this prediction 
are expected . In the Standard Model , a single phase 
in the CKM matrix governs CP violation [3( , but if heavy 
non-SM particles appear in additional penguin diagrams, 
new CP- violating phases could enter and S p o K o would 
not equal sin 2/3 0. Observation of a significant discrep- 
ancy would be a clear signal of new physics. 

In this Letter we present the first observation of the de- 
cay B° — > p°K a and a measurement of the CP-violating 
asymmetries S p o K o and C p o K o from a time-dependent 
maximum likelihood analysis. A non-zero value of SpO^o 
indicates CP violation due to the interference between 
decays with and without mixing. Direct CP violation 
leads to a non-zero value of C p ox° s ■ We take a quasi- two- 
body (Q2B) approach, restricting ourselves to the region 
of the B° — ► 7r + 7r~A"° Dalitz plot dominated by the p° 
and treating other B° — ► tt + it~ K® contributions as non- 
interfering background. The effects of interference with 
other resonances are estimated and taken as systematic 
uncertainties. 

The data were collected with the BABAR detector at the 
PEP-II asymmetric-energy e + e~ storage ring at SLAC. 
An integrated luminosity of 205 fb -1 , corresponding to 
227 x 10 6 BB pairs, was collected at the T(45) reso- 
nance (center-of-mass (CM) energy y/s = 10.56 GeV), 
and 16 fb -1 was collected about 40 MeV below the res- 
onance (off-resonance data). The BABAR detector is de- 
scribed in detail elsewhere 0- Charged particles are de- 
tected and their momenta measured by the combination 
of a silicon vertex tracker (SVT), consisting of five lay- 
ers of double sided detectors, and a 40-layer central drift 
chamber (DCH), both operating in the 1.5 T magnetic 
field of a solenoid. Charged-particle identification is pro- 
vided by the average energy loss in the tracking devices 
and by an internally reflecting ring-imaging Chercnkov 
detector (DIRC) covering the central region. 

We reconstruct B° — ► p°K® candidates (B° ec in the 
following) from combinations of p° and candidates, 
both reconstructed in their tt + tt~ decay mode. For the 
7r + 7r~ pair from the p° candidate, we remove tracks iden- 
tified as very likely to be electrons, kaons, or protons. 



The mass of the p candidate is restricted to the interval 
0.4 < m(7r+7r-) < 0.9GeV/c 2 . The K° s candidate is re- 
quired to have a mass within 13 MeV/c 2 of the nominal 
K® mass [a] and a decay vertex separated from the p° 
decay vertex by at least three times the estimated seper- 
ation measurement uncertainty. In addition, the cosine 
of the angle in the lab frame between the flight di- 
rection and the vector between the p° decay vertex and 
the K® decay vertex must be greater than 0.995. Vetoes 
against B° -> D+ir- and B° -► K*w~(K* -> K°ir+) 
are imposed by requiring that the invariant masses of 
both K®ir combinations are more than 0.055 GeV/c 2 and 
0.040 GeV/c 2 from the K* + and D + masses respec- 
tively. To exclude events with poorly reconstructed ver- 
tices we require the estimated error on At to be less than 
2.5 ps and that |Ai| must be less than 20 ps, where At is 
the proper time difference between the decay of the recon- 
structed B meson (B° ec ) and its unreconstructed partner 
(^tag); ^rec ~ ^tag- It is determined from the measured 
relative displacement of the two B-decay vertices and the 
known boost of the e + e~ system. 

Two kinematic variables are used to discriminate be- 
tween signal and combinatorial background. The first is 
AE, the difference between the measured CM energy of 
the B candidate and y/s/2, where y/s is the CM beam 
energy. The second is the beam-energy substituted mass 
toes = \J (s/2 + p., ■ p B ) 2 /B 2 - p|", where the B r ° cc mo- 
mentum ps and the four- momentum of the initial T(4S) 
state (Bj, pj) are defined in the laboratory frame. We 
require |AB| < 0.15 GeV and 5.23 < m ES < 5.29GeV/c 2 . 

Continuum e + e~ — > qq (q = u,d,s,c) events are the 
dominant background. To enhance discrimination be- 
tween signal and continuum, we use a neural network 
(NN) to combine five variables: the cosine of the angle 
between the B° ec direction and the beam axis in the CM, 
the cosine of the angle between the thrust axis of the 
B° ec candidate and the beam axis, the sum of momenta 
transverse to the direction of flight of the B° ec , and the ze- 
roth and second angular moments Bo. 2 of the energy flow 
about the B° ec thrust axis. The moments are defined by 
Lj = J2iPi x | cos 0i | J ', where Pi is its momentum and 
0i is the angle with respect to the Bj? ec thrust axis of 
the track or neutral cluster i excluding the tracks that 
make up the B° ec candidate. The NN is trained with 
off- resonance data and Monte Carlo (MC)0 simulated 
signal events. 

selected signal events are reconstructed incorrectly 
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with low momentum tracks from the other B meson be- 
ing used to form the p° candidate. In total, 20,073 events 
pass all selection criteria in the on-resonance sample. 

An unbinned extended maximum likelihood fit is used 
to extract the p°K® CP asymmetry and branching frac- 
tion. There are ten components in the fit: signal, contin- 
uum background and eight separate backgrounds from B 
decays. Large samples of MC-simulated events are used 
to identify these specific B backgrounds. Where an in- 
dividual decay mode makes a significant contribution to 
the dataset (one or more events expected in the data) we 
include it as a separate contribution to the fit. Probabil- 
ity density functions (PDFs) are taken from simulation 
with the expected number of B background events fixed 
to values estimated from known branching fractions 
and MC efficiencies (Table HJ. Where only upper limits 
are available, decay modes are not included in the default 
fit but are used in alternate fits to evaluate systematics. 

Events from B decays that do not come from individ- 
ually significant channels are collected together into two 
"bulk" B contributions to the fit (B° and B+). The 
assumption is made that B° — > f o (Q00)K° can be ne- 
glected, with support from 0, which do not require 
this mode to describe B + — > K + tt + it~. 

The events in the data sample have their unrecon- 
structed Bs flavor-tagged as B° or B° with the method 
described in 0. Events are separated into four flavor- 
tagging categories and an "untagged" category, depend- 
ing upon the method used to determine the flavor. Each 
category has a different expected purity and accuracy of 
tagging. The likelihood function for the Nk candidates 
in flavor tagging category k is 



C k = e-^nS {Nse k [(1 - « ul ;Ki: : + f^g 

+Nc,^ k + .\:, :<u 7'» k |. (I) 

where N' k is the sum of the signal and background yields 
for events tagged in category k, N$ is the number of 
p°K® signal events in the sample, is the fraction of 
signal events tagged in category k, is the fraction of 
mis-reconstructed (MR) signal events in tagging category 
k and the superscript CR implies correctly reconstructed 
signal. Nc,k is the number of continuum background 
events that are tagged in category k, and iVeje^k is the 
number of B-background events of class j that are tagged 
in category k. The ^-background event yields are fixed 
in the default fit to values shown in Table [I] The values 
€k and f k are determined from MC for ^-backgrounds 
and from a sample of B decays of known flavor for signal. 
The total likelihood C is the product of the likelihoods 
for each tagging category. 

Each signal and background PDF is defined as: Vu = 
V (m ES ) • V ( AE) ■ V k (NN) • V (cos 9^ ) ■ V ( At ) • V (m I+ w - ) 



Background Mode 


expected 


Bulk B+ 


197±98 


Bulk B" 


197±98 


B° D+7T- 


40±6 


B" — > ri'K° s 




B° -» f (980)K s 
B° -> A- *(1430) + 7r- 


22±4 
7±1 


B° -» p°K* 


3±3 


B° (J&t+Onb 


2±1 



TABLE I: Expected number of events from each B back- 
ground source. 



where toes, AE, NN, m(7r + 7r~) are the variables de- 
scribed previously, and cos 9^+ is the angle between the 
and the 7r + from the p° in the p° meson's center-of- 
mass frame. 

The At PDF for signal events is defined as 



V{At) 



e -|At]/r B 

~ B X 

1 + + li D ) ( S p°ko sin(Am d At)- 
C p o K o cos(Am d At)) ® R Big (At,aAt), (2) 



where tb and Arrid are the average lifetime and eigen- 
state mass difference of the neutral B meson, q = 
+1 (-1) when B° cc = B° (B°), (D) describes the dilu- 
tion effect from imperfect flavor tagging, and AD is the 
difference in this dilution between B° and B° tags. This 
formalism is found to effectively describe both correctly 
and incorrectly reconstructed signal. (D), AD and the 
At resolution function, R a i s (At, <TAt), have parameters 
fixed to values taken from a sample where Bs of known 
flavor can be reconstructed 01 ■ "Untagged" events have 
a (D) of 0, reflecting the lack of tag information. 

The toes, AE, NN, costf^ and to^tt^) PDFs for 
signal and B background are taken from MC simulation. 
In general they are non-parametric, with the exception 
of toes an d AE for signal signal PDFs appear as solid 
curves in Figure^jThe CP parameters for rf ' and foK® 
backgrounds are fixed to C — and S = sin 2/3 (for 
rj' K®) and S = — sin 2/3 (for faK®), in accordance with 
SM expectations. For the remaining B backgrounds the 
parameters C and S are fixed to 0. The PDF parameters 
describing the continuum background are either allowed 
to vary freely in the fit or else determined separately from 
off-resonance data. 

There are 16 free parameters in the fit: the yield of 
signal events, S pK o and C pK o and 13 that parameterize 
the continuum background. The continuum parameters 
are: the yields (5), and those asociated with the second 
order polynomial describing the AE distribution (2), the 
ARGUS [llj function describing the toes distribution (1) 
and the double Gaussian used to model the At distribu- 
tion (5). 

The fit yields 111 ± 19 signal events. We calculate 
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the branching fraction from the measured signal yield, 
efficiency (including the p° — > ir + ir~, K° — > K® and 
Kg — > tt + it~ branching fractions), and the number of 
BB events. The result is B(B° -> p°K°) = (4.9 ± 
0.8 ± 0.9) x 10~ 6 , where the first error is statistical and 
the second systematic. The likelihood ratio between the 
fit result of 111 signal events and the null hypothesis 
of zero signal shows that this is excluded at the 8.7 a 
level. When additive systematic effects are included we 
exclude the null hypothesis at the 5.0a level. The fit for 
CP parameters gives S p o K « = 0.20 ± 0.52 ± 0.24 and 
C a o K o = 0.64 ± 0.41 ±0.20. S 

Figure shows s Vlots 12] of the discriminating vari- 
ables in the fit. Knowledge of the level of background 
and our ability to distinguish it from signal can be gained 
from the errors in these plots. In addition, Fig.HJf) shows 
the ratio £s/(Cs + £b) for all events, where Ls and Lb 
are the likelihoods for each event to be signal or back- 
ground, respectively. 

Figure |21 shows gPlots of At. Untagged events are re- 
moved, and events are split into B® ag tags and B® ag tags. 



)/(N B o+N S o) 



An s Vlot of asymmetry (N B o — N 
as a function of At is also shown 

Systematic errors are listed in Table [H] We estimate 
biases due to the fit procedure from fits to a large num- 
ber of simulated experiments. We vary parameters fixed 
in the nominal fit by their uncertainty and include the 
change in result as the corresponding systematic error. 
The systematic uncertainties arise from sources including 
the parameterization of the signal At resolution function, 
the mistag fractions, and discrepancies between data and 
the simulation including the effect of alternative models 
for resonances. 

We estimate the systematic uncertainly due to ne- 
glecting the intereference between B° — > K ( *tt + it~ from 
both parameterized and full simulations that take in- 
terference into account. We include contibutions from 
p°(770)K°, / (980)ir°, iq(1430)+7r-, K*(892)+ir- and 
/2(1270)id?, as well as two K®irir non-resonant contribu- 
tions. We simulate many samples with different relative 
phases between modes. We also vary the amplitude of 
each mode within limits based on the best available in- 
formation 8, 9J. Each simulation is then subjected to the 
standard selection and fitting procedure. The systematic 
uncertainty is taken from the width of a Gaussian fitted 
to the distribution of the results. 

In summary, we have established the existance of the 
decay B — > p°K° and measured its branching fraction 
with the significance of 5 standard deviations. Our mea- 
surement agrees within errors with B(B° — ► ljK°) as 
measured in |l3| , as expected if a single penguin ampli- 
tude dominates these decays. We have extracted the CP 
violating parameters S and C for B° — » p°K® which are 
consistent with those measured in charmonium channels 





«60 r 
a 



(b) 



340- 



5.24 5.25 5.26 5.27 5.28 5.29 
m ES (GeV/c 2 ) 




0.15 -0.1 -0.05 0.05 0.1 0.15 

AE (GeV) 




0.7 0.8 0.9 

m ra (GeV/c 2 ) 



0.2 0.4 0.6 0.8 1 
L, /(L s . +L„. ) 

Sig Sig Bkg 7 



FIG. 1: sVlots of Maximum Likelihood fit discriminating 
variables: (a) mES, (b) Ai?, (c) Neural Network output, (d) 
cos#,„-+, (e) invariant mass of the 7r + 7r~ combination. Lines 
are projections of signal PDFs for each variable, (f ) is a plot 
of the likelihood of an event being signal calculated for all 
events in our dataset and compared to the predictions of our 
PDF (predicted continuum in light grey, B Background in 
dark grey and signal in white). 
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FIG. 2: s Vlots of At, overlaid with projected signal PDFs, 
split into (a) Stag tags, (b) Stag tags and (c) the asymmetry 
(N B o - N n o )/{N B a + N n o ) as a function of At . 

V "tag ""tag" V "tag ""tag 7 



Mis-reco'd events and fit bias 0.12 0.09 10 

PDF uncertainties 0.13 0.18 2 

Tagging parameters 0.02 0.01 - 

Neglect of intereference 0.14 0.09 7 

p° mass shape 0.07 0.05 3 

B Background BF 0.02 0.10 13 

CP of background 0.04 0.00 - 
Tracking efficiency & B counting - - 6 



* Also at Laboratoire de Physique Corpusculaire, 
Clermont-Ferrand, France 

t Also with Universita di Perugia, Dipartimento di Fisica, 
Perugia, Italy 

* Also with Universita della Basilicata, Potenza, Italy 

[1] BABAR Collaboration, B. Aubert et al, Phys. Rev. Lett. 

94, 161803 (2005), BELLE Collaboration, K. Abe et al, 

Phys. Rev. D 71, 072003 (2005). 
[2] G. Buchalla et al, JHEP 0509, 074 (2005); M. Beneke, 

Phys. Lett. B 620, 143 (2005); H. Y. Cheng, C. K. Chua 

and A. Soni, Phys. Rev. D 72, 014006 (2005). 
[3] N. Cabibbo, Phys. Rev. Lett. 10, 531 (1963); 

M. Kobayashi and T. Maskawa, Prog. Theor. Phys. 49, 

652 (1973). 

[4] Y. Grossman and M. P. Worah, Phys. Lett. B 395, 241 
(1997); M. Ciuchini et al, Phys. Rev. Lett. 79, 978 
(1997); D. London and A. Soni, Phys. Lett. B 407, 61 
(1997). 

[5] BABAR Collaboration, B. Aubert et al, Nucl. Instrum. 
Methods Phys. Res., Sect. A A479, 1 (2002). 

[6] Particle Data Group, S. Eidelman et al, Phys. 
Lett. B 592, 1 (2004). 

[7] S. Agostinelli et al. [GEANT4 Collaboration], Nucl. In- 
strum. Meth. A 506, 250 (2003). 

[8] BABAR Collaboration, B. Aubert et al, Phys. Rev. D 72 
072003 (2005). 

[9] BELLE Collaboration, A. Garmash et al, Phys. Rev. D 
71, 092003 (2005). 
[10] BABAR Collaboration, B. Aubert et al, Phys. Rev. D 66, 
032003 (2002). 

[11] ARGUS Collaboration, H. Albrecht et al, Phys. Lett. 

B241, 278 (1990). 
[12] M. Pivk and F. R. Le Diberder, Nucl. Instrum. Methods 

Phys. Res. Sect. A A555, 356 (2005). 
[13] BABAR Collaboration, B. Aubert et al, Phys. Rev. D 74, 

011106 (2006). 



Total 



0.24 0.20 19 



TABLE II: Summary of contributions to the systematic error. 



